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Abstract 

We review the present status for the determinations of the light and heavy quark masses, the light 
quark chiral condensate and the decay constants of light and heavy-light (pseudo)scalar mesons from 
QCD spectral sum rules (QSSR). Bounds on the light quark running masses at 2 GeV are found to 
be (see Tables @ and |): 6 MeV< (m d +fn u )(2) <11 MeV and 71 MeV< m s (2) <148 MeV. The 
agreement of the ratio m s /(m u + md) = 24.2 in Eq. (IB"?] ) from pseudoscalar sum rules with the 
one (24.4 ± 1.5) from ChPT indicates the consistency of the pseudoscalar sum rule approach. QSSR 



predictions from different channels for the light quark running masses lead to (see Section 10.3 ): 
m s (2) = (117.4±23.4) MeV, (m d +m u )(2) = (10.1±1.8) MeV, (m d -W u )(2) = (2.8±0.6) MeV with the 
corresponding values of the RG invariant masses. The different QSSR predictions for the heavy quark 
masses lead to the running masss values: rn c (m c ) = (1.23±0.05) GeV and m^mj,) = (4.24±0.06) GeV 
(see Tables || and ^), from which one can extract the scale independent ratio mi,/m s = 48.8 ± 9.8. 
Runned until Mz, the 6-quark mass becomes: m^Mz) = (2.83 ± 0.04) GeV in good agreement with 
the average of direct measurements (2.82 ± 0.63) GeV from three-jet heavy quark production at LEP, 
and then supports the QCD running predictions based on the renormalization group equation. As a 
result, we have updated our old predictions of the weak decay constants /v(i.3)j /k'(i.46); /a (o.98) 
and /if*(i.43) (see Eqs. ( |S7| ) and (|S9"|)). We obtain from a global fit of the light (pseudo)scalar and B s 
mesons, the flavour breakings of the normal ordered chiral condensate ratio: (ss)/(uu) = 0.66 ± 0.10 



(see Eq. (113)). The last section is dedicated to the QSSR determinations of /d (s1 and Jb, 
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lr This review updates and completes the reviews jjj, E| and some parts of the book jjj. It has been extracted from a 



chapter of the forthcoming book: QCD as a theory of hadrons: from partons to confinement 



1 Introduction 



One of the most important parameters of the standard model and chiral symmetry is the light and 
heavy quark masses. Light quark masses and chiral condensates are useful for a much better under- 
standing of the realizations of chiral symmetry breaking and for some eventual explanation of 
the origin of quark masses in unified models of interactions TH . Within some popular parametrizations 
of the hadronic matrix elements || , the strange quark mass can also largely influence the Standard 
Model prediction of the CP violating parameters e'/e, which have been mesured recently How- 
ever, contrary to the QED case where lcptons are observed, and then the physical masses can be 
identified with the pole of the propagator (on-shell mass valuc)[| the quark masses are difficult to 
define because of confinement which docs not allow to observe free quarks. However, despite this 
difficulty, one can consistcnly treat the quark masses in perturbation theory like the QCD coupling 
constant. They obey a differential equation, where its boundary condition can be identified with the 
rcnormalized mass of the QCD lagrangian. The corresponding solution is the running mass, which is 
gauge invariant but renormalization scheme and scale dependent, and the associated renrmalization 
group invariant mass. To our knowledge, these notions have been introduced for the first time in 
[Xj]. In practice, these masses are conveniently defined within the standard MS-scheme discussed in 
previous chapters. In addition to the determination of the ratios of light quark masses (which are 
scale independent) from current algebra and from chiral perturbation theory (ChPT), its modern 
version jl{|-|l9), a lot of effort reflected in the literature [^0| has been put into extracting directly 
from the data the running quark masses using the SVZ |2l[] QCD spectral sum rules (QSSR) j|, LEP 
experiments and lattice simulations. The content of these notes is: 

• a review of the light and heavy quark mass determinations from the different QCD approaches. 

• a review of the direct determinations of the quark vacuum condensate using QSSR and an update 
of the analysis of its flavour breakings using a global fit of the meson systems. 

• An update of the determinations of the light (pseudo)scalar decay constants, which, in particular, 
are useful for understanding the qq contents of the light scalar mesons. 

• A review of the determinations of the weak leptonic decay of the heavy- light pseudoscalar mesons 
D( s) and S (s) . 

This review develops and updates the review papers [Q, ||] and some parts of the book [||. It also 
updates previous results from original works. 



2 Definitions of perturbative quark masses in QCD 

Let's remind the meaning of quark masses in QCD. One starts from the mass term of the QCD 
lagrangian: 

C m = triilpilpi , (1) 

where mi and ?pi are respectively the quark mass and field. It is convenient to introduce the dimen- 
sionless mass Xi{v) = mi(v)/is, where v is the renormalization scheme subtraction constant. The 
running quark mass is a solution of the differential equation: 

= (1 + 7 (a.))z<(t) : Xiit = 0) = Xi(v) . (2) 

2 For a first explicit definition of the perturbative quark pole mass in the MS-scheme, sec [jll], [l2j (renormalization- 
scheme invariance) and [113] (regularization-scheme invariance). 
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In the MS-scheme, its solution to order a 3 (a s = a s /7r) is: 



mi(u) = m 4 (-/3i« s H)" 7l//3l |i + ^(^--|)« s M 
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where 7.; and /?,; arc the 0{a\) coefficients of the quark-mass anomalous dimension and /^-function, 
which read for SU(3) C x SU{n) f : 
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The invariant mass rhi has been introduced for the first time by jlij in connection with the analysis 
of the breaking of the Weinberg sum rules by the quark mass terms in QCD. For the heavy quarks, 
one can also define a perturbative (short distance) pole mass at the pole of the propagator. The IR 
finiteness of the result to order a 2 , has been explicitly shown in jll], 12 1. The independence of M po i e on 
the choice of the regularization-scheme has been demonstrated in [13]. The extension of the previous 



result to order a s is: 22 



Mpole = m(p 2 ) 



Kc 



- + In — ) ( — 

3 m 2 



221 

Ya 



13 

7 

36 



n I In 



15 



n 

12 



2 V 



where in the RHS m is the running mass evaluated at p 2 and: 
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Kq = 17.1514 - 1.04137n +-^A 
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For < r < 1, A(r) can be approximated, within an accuracy of 1% by: 

7T 2 

A(r) ~ — r - 0.597r 2 + 0.230r 3 . 



(5) 
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It has been argued that the pole masses can be affected by nonperturbative terms induced by the 
resummation of the QCD perturbative series |2^] and alternative definitions free from such ambiguities 
have been proposed (residual mass |24| and IS mass ^5|). Assuming that the QCD potential has no 
linear power corrections, the residual or potential-subtracted (PS) mass is related to the pole mass as: 



M PS = M, 



pole 
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(8) 



The IS mass is defined as half of the perturbative component to the 3 Si QQ ground state, which is 
half of its static energy (2M po i e + V) 0. The running and short distance pole mass defined at a given 
order of PT series will be used in the following discussions. 

3 These de finitions might still be affected by a dimcnsion-2 term advocated in | p^ , p?| pj| |, which might limit their 
accuracy |2i3]. 
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3 Ratios of light quark masses from ChPT 



The ratios of light quark masses are well-determined from current algebra || , and ChPT Jl5| . In this 
approach, the meson masses are expressed using a systematic expansion in terms of the light quark 
masses: 



Ml + = {m u + m d )B + <D(m 2 ) + 
(m u + m s )B + 0(m 2 ) + . 

)B + 0(m 2 ) + . 



M 2 K0 



(9) 



where B = —(rpift)/^ from the Gell-Mann, Oakes, Renner relation 



mif 2 ~ -(m u + m d )(#) + 0{mf) 



(10) 



However, only the ratio, which is scale independent can be well determined. To leading order in m 
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Including the next order + electromagnetic corrections, the ratios of masses are constrained on the 
ellipse: 
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where: 
average 
formula 
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m%)(l + 
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'')/(m 2 — m 2 ) — 22.7 ± 0.8 using the value of the r\ 



(12) 



from the PDG 



though this value can well be in the range 22-26, to be compared with the Dashen's 
of 24.2; rh = (l/2)(m u + m d ). In the figure of one shows the range spanned by 
— m)/(m d — m u ) and the corrections to the GMO mass formula Am : Af| = (1/3)(4M|- — 
Am)- The Weinberg mass ratio (j| is also shown which corresponds to the Dashen's formula 
43. At the intersection of different ranges, one deduces 0: 



— = 0.553 ±0.043 , — 

m d m d 



18.9 ±0.8, 



2m s 



(m d + m u ) 



24.4 ± 1.5. 



(13) 



The possibility to have a m u = advocated in |33|] appears to be unlikely as it implies too strong 
flavour symmetry breaking and is not supported by the QSSR results from 2-point correlators of the 
divergences of the axial and vector currents, as will be shown in the next sections. 



4 Bounds on the light quark masses 

In QSSR, the estimate and lower bounds of the sum of the light quark masses from the pseudoscalar 
sum rule have been firstly done in |34|, |35| , while a bound on the quark mass difference has been firstly 
derived in |3(|. The literature in this subject of light quark masses increases with time ^] However, 
it is in some sense quite disappointing that in most of the published papers no noticeable progress 

4 In Generalized ChPT, the contribution of the m 2 -term can be as large as the m one Jj]||, which modifies drastically 
these ratios. 

5 This figure will be shown in the book version of the paper. 
6 Previous works are are reviewed in P, bl- 
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has been done since the previous pioneering works. The most impressive progress comes from the 
QCD side of the sum rules where new calculations have become available both on the perturbative 
radiative corrections known to order o? s 0, [57| [38| and on the nonperturbative corrections [plf-|4l| 
fl. Another new contribution is due to the inclusion of the tachyonic gluon mass as a manifestation 
of the resummation of pQCD series |26|, [2?], ^8|. Alas, no sharp result is available on the exact size 
of direct instanton contributions advocated to be important in this channel fti| , while [fill claims 
the opposite. Though the instanton situation remains controversial, recent analysis |44|, fusing the 
results of |46) based on the ILM of jfTj indicates that this effect is negligible justifying the neglect of 
this effect in different analysis of this channel. However, it might happen that adding together the 
effect of the tachyonic gluon to the one of direct instanton might also lead to a double counting in a 
sense that they can be two alternative ways for parametrizing the nonperturbative vacuum p8| . In 
absence of precise control of the origin and size of these effects, we shall consider them as new sources 
of errors in the sum rule analysis. 



4.1 Bounds on the sum of light quark masses from pseudoscalar channels 

Lower bounds for (jn u + rrid) based on moments inequalities and the positivity of the spectral functions 
have been obtained, for the first time, in |54| [35| . These bounds have been rederived recently in |l^, 
to order a s . As cheked in jjj for the lowest moment and redone in [ [i"5"| for higher moments, the inclusion 
of the Qg term decrease by about 10 to 15% the strength of these bounds, which is within the expected 
accuracy of the result. 

For definiteness, we shall discuss in details the pseudoscalar two-point function in the us channel. The 
analysis in the ud channel is equivalent. It is convenient to start from the second derivative of the 
two-point function which is superficially convergent: 

f°° 9 1 

*"W 2 )=/ dt U I ^3 ~ Im *ffl- ( 14 ) 







(t + Q 2 



The bounds follow from the restriction of the sum over all possible hadronic states which can contribute 
to the spectral function to the state(s) with the lowest invariant mass. The lowest hadronic state which 
contributes to the corresponding spectral function is the K-pole. From eq. ([H]) we then have 

2 f°° 2 1 

*' 5 ' (Q2) = WFTW 2/iA4 + L d W^ Ms(t) ' (15) 

where t = (Mk + 2m w ) 2 is the threshold of the hadronic continuum. 

It is convenient to introduce the moments £jv(<3 2 ) of the hadronic continuum integral 

One is then confronted with a typical moment problem (see e.g. Ref. pOj.) The positivity of the 
continuum spectral function ^Im^s^) constrains the moments £at(Q 2 ) and hence the l.h.s. of Eq. 



(15) where the light quark masses appear. The most general constraints among the first three moments 
for N = 0,1,2 are: 

E (Q 2 )>0, S!(Q 2 )>0, S 2 (Q 2 )>0; (17) 

S (Q 2 )-S!(Q 2 )>0, £i(Q 2 )-£ 2 (Q 2 )>0; (18) 

S (Q 2 )S 2 (Q 2 )- (Si(Q 2 )) 2 >0. (19) 

The inequalities in Eq. ( |l^ ) are in fact trivial unless 2Q 2 < to, which constrains the region in Q 2 to 
too small values for pQCD to be applicable. The other inequalities lead however to interesting bounds 

7 See also the chapter on two-point functions where more references to original works are given. 



4 



which we next discuss. 

The inequality Eo(<3 2 ) > results in a first bound on the running masses: 



[m s (Q 2 ) + m u {Q 2 )]'' 



> 



16ir 2 2f%M* 



1 



Q 4 



(20) 



where the dots represent higher order terms which have been calculated up to O(o^), as well as 
non-perturbative power corrections of O (l/Q 4 ) and strange quark mass corrections of O (m 2 /Q 2 ) 
and O (m 4 /<2 4 ) including 0(a s ) terms . Notice that this bound is non-trivial in the large-iV c limit 



{fx ~ C(-^c)) and in the chiral limit (m s ~ Mk)- The bound is of course a function of the choice of 
the euclidean Q-value at which the r.h.s. in Eq. ( |20| ) is evaluated. For the bound to be meaningful, 
the choice of Q has to be made sufficiently large. In Ref. jl8) it is shown that Q > 1.4 GeV is already 
a safe choice to trust the pQCD corrections as such. The lower bound which follows from Eq. fl20|) for 
m u + m s at a renormalization scale /i 2 = 4 GeV 2 results in the solid curves shown in Fig. |l| below. 



120n 



^100 




1.8 2.0 
Q [GeV] 



Figure 1: Lower bound in MeV to order a s for (m s + m„)(2) versus Q in GeV from Eq. ( |20| ) for 
A 3 = 290 MeV (upper curve) and 380 MeV (lower curve). Quark mass values below the solid curves 
in Fig. |l| are forbidden by the bounds. 



The resulting value of the bound at Q = 1.4 GeV is: 

(m s + m u )(2) > 80 MeV (m u + m d )(2) > 6.6 MeV , 



(21) 



if one uses either ChPT and the previous SR analysis for the mass ratios. Radiative corrections tend 
to decrease the strengths of these bounds. Their contributions to the second moment of the two-point 
function are (see previous part of the book) : 



3 (m u + m s y 
8^ Q 2 



1 + I — ) -M IT*) ( — 



77.368 I — 



At this scale, the PT series converges quite well and behaves as: 

Parton 1 + 0.45 + 0.22 + 0.15 . 

Including these higher order corrections, the bounds become: 

(m s + m u )(2) > (71.4 ± 3.7) MeV => (m u + m d )(2) > (5.9 ± 0.3) MeV , 



(22) 

(23) 
(24) 
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The bound will be saturated in the extreme limit where the continuum contribution to the spectral 
function is neglected. 

The quadratic inequality in Eq. (|l^) results in improved lower bounds for the quark masses which we 
show in Fig. || below. 




3.0 
Q [GeV] 



Figure 2: The same as in Fig. [I] but from the quadratic inequality to order a s . 



The quadratic bound is saturated for a <5-like spectral function representation of the hadronic contin- 
uum of states at an arbitrary position and with an arbitrary weight. This is certainly less restrictive 
than the extreme limit with the full hadronic continuum neglected, and it is therefore not surprising 
that the quadratic bound happens to be better than the ones from Eat(Q 2 ) for N — 0, 1, and 2. Notice 
however that the quadratic bound in Fig. ||| is plotted at higher Q-values than the bound in Fig. [3]. 
This is due to the fact that the coefficients of the perturbative series in a s (Q 2 ) become larger for the 
higher moments. In Ref Es| it is shown that for the evaluation of the quadratic bound Q > 2 GeV is 
already a safe choice. 

Similar bounds can be obtained for (m u + m<j) when one considers the two-point function associated 
with the divergence of the axial current 



d (1 A ,i (x) = (m d + m„) :d(x)ij 5 u(x): . 



(25) 



The method to derive the bounds is exactly the same as the one discussed above and therefore we 
only show, in Fig. [5] below, the results for the corresponding lower bounds which one obtains from the 
quadratic inequality. At Q — 2 GeV, one can deduce the lower bounds from the quadratic inequality: 

(m s + m„)(2) > 105 MeV , (m u + m d ){2) > 7 MeV (26) 

The convergence of the QCD series is less good here than in the lowest moment. It behaves as p5[: 
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which numerically reads: 



3 \ 7T 
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1 + 0.83 + 0.61 + 0.51 . 
This leads to the radiatively corrected lower bound to order a^: 

(m s + m„)(2) > (82.7 ± 13.3) MeV , (m„ + m d )(2) > (6 ± 1) MeV 



(27) 
(28) 
(29) 
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Figure 3: Lower bound in MeV for (rrid + m u )(2) from the quadratic inequality to order a s 



where the error is induced by the truncation of the QCD series which we have estimated to be about 
the contribution of the last known a 3 term of the series 0. From the previous analysis, and taking into 
account the uncertainties induced by the higher order QCD corrections, the best lower bound comes 
from the lowest inequality and is given in Eq. ( ^4|) . The result is summarized in Table ^. 

Table 1: Lower bounds on m Ui( 2 jS (2) in MeV 
Observables Sources Authors 



rn u + m d 
6 

6.8 



7T LRT97J1, Y97gg] (updated here to order a 3 s ) 

(V'V)+GMOR DN98(5j (leading order) 



md - m u 
1.1 



Kir 



Y97W3 (updated here to order a 3 s 



71.4 

90 



K LRT97_B| (updated here to order a 3 s ) 

(V>V)+ChPT DN98(5| (leading order) 



8 In Ref. |4^, alternative bound has been derived using a Holder type inequality. The lower bound obtained from 
this method, which is about 4.2 MeV is weaker than the one obtained previously. 
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4.2 Lower bound on the light quark mass-difference from the scalar sum 
rule 

As in p6| , one can extract lower bound on the light quark mass difference (m u — m<j) and (m u — m s ) 
working with the two-point function associated to the divergence of the vector current: 

d n v u q = ( m « _ m q) '■ $u(i)i>q ■ ■ (30) 

The most recent analysis has been done in Q . We have updated the result by including the ct^-term. 
It is given in Table fjj 



Table 2: Upper bounds on m u> d tS (2) in MeV. 



Observables Sources Authors 



m„ + nid 

11.4 ( ; 0V)+GMOR DN98|| (leading order) 



148 (?/>V)+ChPT DN98|p| (leading order) 

147 ± 21 e+e~ + r-decay SN99|52| (to order a 3 s ) 



4.3 Bounds on the sum of light quark masses from the quark condensate 
and e + e — > / = hadrons data. 

Among the different results in |5l| , we shall use the range of the chiral (ijnp) = (uu) ~ (dd) condensate 
from the vector form factor of D — > K*lv. Using three-point function sum rules, the form factor reads 
to leading order: 

jr/A\ m c (m D + m K «) u 2 2^ , 2 1 

7(0) = — — g exp[(m D -m c )n + m K ,T 2 \ (31) 

x W)(-l + Mo 2 (-^ + ^ + 2 ^-^ TO V 1 r 2 ) 



7? 



16tt - 2m c 
._a.p<^)(— nr,-— (1 

2m"? — m 2 m s 9 —m c 2 2m s ? 2 4 T2 

36 9 9 9 9 m r 



e m c Tl PS20 
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ds 2 [ 10 d Sl p v ( Sl , S 2)e- SlTl - S2T2 

J S2 -\-m^ 



with Pv(si,S 2 ) = , 9 , 3 rrr X (32) 

47T (Si - S2) 13 

|m s ((si + s 2 )(si - m^) - 2sis 2 ) + m c ((s 1 + s 2 )s 2 - 2s 2 (si - m 2 c ))\ . 

The factor p ~ 2 ~ 3 expresses the uncertainty in the factorization of the four quark condensate. 
In our numerical analysis, we start from standard values of the QCD parameters and use fx* — 
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0.15 GeV(/ 7r = 93.3 MeV). The value of fjj ~ (1.35 ± 0.07) /„. is consistently determined by a two- 
point function sum rule including radiative corrections as we shall see in the next chapter, where 
the sum rule expression can, e.g., be found in ||. The following parameters enter only marginally: 
m s (l GeV) = (0.15 - 0.19) GeV, s 10 = (5 ~ 7) GeV 2 , s 20 = (1.5 ~ 2) GeV 2 . Using the conservative 
range of the charm quark mass: m c (pole) between 1.29 and 1.55 GeV (the lower limit comes from the 
estimate in || and the upper limit is 1/2 of the J/ty mass), one can deduce the running condensate 
value at 1 GeV |lj: 

0.6 < (i/iip)/[-225 MeV] 3 < 1.5. (33) 

This result has been confirmed by the lattice [[53| . Using the GMOR relation: 

2m 2 / 2 = -{m u + m d )(uu + dd) + C(m 2 ) . (34) 

one can translate the upper bound into a lower bound on the sum of light quark masses. The lower 
bound on the chiral condensate can be used in conjonction with the positivity of the m 2 correction in 
order to give an upper bound to the quark mass value. In this way, one obtains: 

6.8 MeV < (m„ + m d )(2 GeV) < 11.4 MeV . (35) 

The resulting values are quoted in Tables [l] and |^. We expect that these bounds are satisfied within 
the typical 10% accuracy of the sum rule approach. 

We also show in Table || the upper bound obtained in []52| by using the positivity of the spectral 
function from the analysis of the e + e~ — > I = hadrons data where the determination will be 
discussed in the next section. 



5 Sum of light quark masses from pseudoscalar sum rules 
5.1 The (pseudo) scalar Laplace sum rules 

The Laplace sum rule for the (pseudo) scalar two-point correlator reads (see e.g. Q-[Ell: 



/' 

Jo 



dt exp (—tr) — Imty^-j (t) ~ (m u ± m^) 2 



8tt 2 



(i-Pi)(i+4 0) ) + £4 n) 



n=2 



(36) 



where the indices 5 and + refer to the pseudoscalar current. Here, r is the Laplace sum rule variable, 
t c is the QCD continuum threshold and mi is the running mass to three loops, 

p x = (l + t c r)exp(-t c r) . (37) 

Using the results compiled in previous Chapter , the perturbative QCD corrections read for n flavours: 

11 
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where CV=4/3 and je — 0.5772... is the Euler constant; 71, 72 and ft, ft are respectively the mass 
anomalous dimensions and ftfunction coefficients defined in a previous chapter. For three colours and 
three flavours, they read: 



71 = 2 , 72 = 91/12 , ft = -9/2 , fa 
In practice, the perturbative correction to the sum rule simplifies as: 



r(0) 



4.82a s + 21.98af + 53.14^ + O(a^) 



(39) 
(40) 



Introducing the RGI condensates defined in the previous chapter , the non perturbative contributions 
are lEl: 
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— + C F — + - 7iJ + - (4 7l + ft) 
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8tt 2 , 



16 
27^ 



(41) 



Beyond the SVZ expansion, one can have two contributions: 
• The direct instanton contribution can be obtained from 

~2 r 



and reads: 



ynst = £± exp ( _ rc) y K ^ c) + Ri{rc) 



(42) 



with: 



Pc/(2t); p c w 1/600 MeV 1 being the instanton radius; Ki is the Mac Donald 



function. However, one should notice that analogous contribution in the scalar channel leads to 
some contradictions ( |p7| and private communication from Valya Zakharov). 



The tachyon gluon mass contribution can be deduced from [27 



fitach _ _ 4 



A 



(43) 



where (a s /7r)A 2 ~ -0.06 GeV 2 ||, 
which completes the different QCD contributions to the two-point correlator. 
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5.2 The ud channel 



From the experimental side, we do not still have a complete measurement of the pseudoscalar spectral 
function. In the past Q, one has introduced the radial excitation tt' of the pion using a NWA where 
the decay constant has been fixed from chiral symmetry argument |Q and from the pseudoscalar sum 
rule analysis itself ]56|,||5^|,||3[, through the quantity: 



(44) 



Below the QCD continuum t c 
radial excitation and reads: 



1 



the spectral function is usually saturated by the pion pole and its first 



dtexp(-tT)-Im* s (i) ~ 2mlf* exp (-m 2 r) 



1 + r n exp 



Mi 



(45) 



The theoretical estimate of the spectral function enters through the not yet measured ratio r v 
discussions of the sum rule analysis can be found in |5 



55, 5 



Detailed 

However, this channel is quite peculiar 
due to the Goldstone nature of the pion, where the value of the sum rule scale (1/r for Laplace and 
t c for FESR) is relatively large of about 2 GeV 2 compared with the pion mass, where the duality 
between QCD and the pion is lost. Hopefully, this paradox can be cured by the presence of the new 
1/q 2 [^6[ H|, |28| due to the tachyonic gluon mass, which enlarges the duality region to lower scale and 
then minimizes the role of the higher states into the sum rule. This naive NWA parametrization has 
been improved in |5§| l by the introduction of threshold effect and finite width corrections. Within the 
advent of ChPT, one has been able to improve the previous parametrization by imposing constraints 
consistent with the chiral symmetry of QCD fl59f. In this way, the spectral function reads: 



1. 



Imtf 5 (f) ~ Imlfi 5(t - mi) + 6(t - 9<) 



1 



t 



(16^ 2 / 2 ) 2 18 



(46) 



with: 




2{t-mlf 



-\{t,u,mj 1 

3 \ u 



10ml 



(t ml) 



[3((u~ml)-t + 10ml] , (47) 



where A (a, b, c) 



2ab — 26c — 2ca is the usual phase space factor. Based on this 



parametrization but including finite width corrections, a recent re-analysis of this sum rule has been 
given to order a 2 |5jJ. Result from the LSR is, in general, expected to be more reliable than the one 
from the FESR due to the presence of the exponential factor which suppresses the high-energy tail 
of the spectral function, though the two analysis are complementary. In |5l| , FESR has been used 
for matching by duality the phcnomenological and theoretical parts of the sum rule. This matching 
has been achieved in the energy region around 2 GeV 2 , where the optimal value of m u + m& has been 
extracted. In the LSR analysis has been updated by including the c? s correction obtained in |38| . 
In this way, we get: 

(m„ + m d ){2 GeV) = (9.3 ± 1.8) MeV , (48) 

where we have converted the original result obtained at the traditional 1 GeV to the lattice choice of 
scale of 2 GeV through: 

mi(l GeV) ~ (1.38 ± 0.06) 7^(2 GeV), (49) 

for running, to order q^, the results from 1 to 2 GeV. This number corresponds to the average value 
of the QCD scale A 3 ~ (375 ± 50) MeV from PDG |(J and f(J. Analogous value of (9.8 ± 1.9) MeV 
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for the quark mass has also been obtained in |6l[] to order as an update of the |39| result. We take 
as a final result the average from § and ||: 

(m„ + m d )(2 GeV) = (9.6 ± 1.8) MeV , (50) 

The inclusion of the tachyonic gluon mass term reduces this value to [ p7| : 

A tach (m u + m d )(2 GeV) ~ -0.5 MeV . (51) 

As already mentioned, adding to this effect the one of direct instanton might lead to a double counting 
in a sense that they can be alternative ways for parametrizing the nonperturbative QCD vacuum. 
Considering this contribution as another source of errors, it gives: 

A ms * (m„ + m d ) (2 GeV) ~ -0.5 MeV . (52) 

Therefore, adding different sources of errors, we deduce from the analysis: 

{m u +m d ){2 GeV) = (9.6 ± 1.8 ±0.4 ±0.5 ±0.5) MeV , (53) 

leading to the conservative result for the sum of light quark masses: 

(M u + m d ){2 GeV) = (9.6 ± 2.0) MeV (54) 

The first error comes from the SR analysis, the second one comes from the running mass evolution 
and the two last errors come respectively from the (eventual) tachyonic gluon and direct instanton 



contributions. This result is in agreement with previous determinations ||, ]55|-[£8|], 62 f33| though 
we expect that the errors given there have been underestimated. One can understand that the new 
result is lower than the old result || |5f| obtained without the a 2 s and a z s terms as both corrections 
enter with a positive sign in the LSR analysis. However, it is easy to check that the QCD perturbative 
series converge quite well in the region where the optimal result from LSR is obained. Combining the 
previous value in Eq. (|J) with the ChPT mass ratio, one can also deduce: 

m s (2 GeV) = (117.1 ± 25.4) MeV . (55) 

5.3 The us channel and QSSR prediction for the ratio m s /(m u + md) 

Doing analogous analysis for the kaon channel, one can also derive the value of the sum [m u + m s ). 
The results obtained from |^5| updated to order oP s and from |64| are shown in Table || given in PJ 
but updated. We add to the original errors the one from the tachyonic gluon (5.5%), from the direct 
instanton (5.5%) and the one due to the evolution from 1 to 2 GeV (4.4%), which altogether increases 
the original errors by 8.9%. Therefore, we deduce the (arithmetic) average from the kaon channel: 

m s (2 GeV) = (116.0 ± 18.1) MeV , (56) 

One should notice here that, unlike the case of the pion, the result is less sensitive to the contribution 
of the higher states continuum due to the relatively higher value of Mk, though the parametrization 
of the spectral function still gives larger errors than the QCD series. It is interesting to deduce from 
Eqs. ( p4[ ) and (|56|), the sum rule prediction for the scale invariant quark mass ratios: 

rs =_?^£_~24.2 > (57) 

where we expect that the ratio is more precise than the absolute values due to the cancellation of the 
systematics of the SR method. This ratio compares quite well with the ChPT ratio |?J: 

rf A = 24.4 ±1.5, (58) 
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and confirms the self-consistency of the pseudoscalar SR approach. This is a non-trivial test of the SR 
method used in this channel and may confirm a posteriori the neglect of less controlled contributions 
like e.g. direct instantons. 



Table 3: QSSR determinations of m s (2) in MeV to order a 3 . Some older results have been updated 
by the inclusion of the higher order terms. The error contains the evolution from 1 to 2 GeV. In 
addition, the errors in the (pseudo)scalar channels contain the ones due to the small size instanton 
and tachyonic gluon mass. Their quadratic sum increases the original errors by 8.9%. The estimated 
error in the average comes from an arithmetic average of the different errors. 



Channels 



m s (2) 



Comments 



Authors 



Pion SR + ChPT 
(#)+ ChPT 



117.1 ±25.4 

129.3 ±23.2 
117.1 ±49.0 



0{al) 



SN99 § Eq. @ 



N, B- B* (l.o) DN98 
D K*lv (l.o) DN98 



Eq. (64) 
Eq. (65) 



Kaon SR 



Scalar SR 



119.6 ± 18.4 updated to 0{p? s ) SN89 



112.3 ± 23.2 0{al) 
116 ±12.8 

148.9 ±19.2 0{al) 
103.6 ± 15.4 
115.9 ±24.0 
115.2 ± 13.0 
99 ± 18.3 



r-like 6 SR: 



-r-decay 129.2 ± 25.6 average 



AS" = — 1 part of r-decay 



-+46.7 
' -57 



169.5: 
144.9 ±38.4 
114±23 
125.7 ±25.4 
115 ± 21 



O 



«2) 




DPS99 
KM01 

CPS97 _ 
CFNP97 j|| 
J98 pT 
M99 
JOP01 



<D(a 3 s ) SN99 




ALEPH99 
CKP98 [|67| 
PP99 J& 
KKP00 
KM01 |0| 
CDGHKK01171 



Average 



117.4 ±23.4 



* Not included in the average. 



6 Direct extraction of the chiral condensate (uu) 



As mentioned in previous section, the chiral uu condensate can be extracted directly from the nucleon, 
B*-B splitting and vector form factor of D — * K*lv, which are particularly sensitive to it and to the 
mixed condensate (-tpa^ ' {\ a /2)G a ^ v \j)) = Mq^i/j) p^ . We have already used the result from the 
D — > K*lv form factor in order to derive upper and lower bounds on (m u + rrid). Here, we shall 
use the informations from the nucleon and from the B*-B splitting in order to give a more accurate 
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estimate. In the nucleoli sum rules [fT2[-||7q], j|, which seem, at first sight, a very good place for 
determining (jpi/i) , we have two form factors for which spectral sum rules can be constructed, namely 
the form factor F\ which is proportional to the Dirac matrix jp and F% which is proportional to 
the unit matrix. In F\ the four quark condensates play an important role, but these are not chiral 
symmetry breaking and are related to the condensate J^ip) only by the factorization hypothesis |2l]] 
which is known to be violated by a factor 2-3 [j72[ ||. The form factor F2 is dominated by 
the condensate (tpip) and the mixed condensate (ipaGip), such that the baryon mass is essentially 
determined by the ratio Mq of the two condensates: 

M 2 = (^aGV)/W) (59) 



Therefore from nucleon sum rules one gets a rather reliable determination of Mq [78 



Mq = (.8 ± .1) GeV 2 . (60) 

A sum rule based on the ratio F2/F1 would in principle be ideally suited for a determination of 
(ipip) but this sum rule is completely unstable |74| due to fact that odd parity baryonic excitations 
contribute with different signs to the spectral functions of F\ and F2. In the correlators of heavy 
mesons (B, B* and D, D*) the chiral condensate gives a significant direct contribution in contrast to 
the light meson sum rules H, since, here, it is multiplied by the heavy quark mass. However, the 
dominant contribution to the meson mass comes from the heavy quark mass and therefore a change 
of a factor two in the value of (^ip) leads only to a negligible shift of the mass. However, from the 
B-B* splitting one gets a precise determination of the mixed condensate (i^aGip) with the value |8(i[| 

{tpaGip) = -(9 ± 1) x 10~ 3 GeV 5 , (61) 

which combined with the value of Mq given in Eq. ( |60|) gives our first result for the value of (ipip) at 
the nucleon scale: 

(ijijj)(M N ) = -[(225±9±9) MeV] 3 , (62) 

where the last error is our estimate of the systematics and higher order contributions. Using the 
GMOR relation, one can translate the previous result into a prediction on the sum of light quark 
masses. The resulting value is: ]5l| : 

(m u + m d ){2 GeV) = (10.6 ± 1.8 ± 0.5) MeV , (63) 

where we have added the second error due to the quark mass evolution. Combining this value with 
the ChPT mass ratio, one obtains: 

m s (2 GeV) ~ 129.3 ± 23.2 MeV . (64) 

Alternatively, one can use the central value of the range given in Eq. ( |33"1) in order to deduce the 
estimate: 

(m u + m d ){2 GeV) = (9.6 ± 4 ± 0.4) MeV =*> m s (2 GeV) ~ (117.1 ± 49.0) MeV . (65) 
The results for m s are shown in Table [|. 

7 Final estimate of (m u + m^) from QSSR and consequences 

on m ui rrid and m s 

One can also notice the impressive agreement of the previous results from pseudoscalar and from the 
other channels. As the two results in Eqs. (|54|), ( |63| ) and ([35]) come from completely independent 
analysis, we can take their geometric average and deduce the final value from QSSR: 

(m u + m d )(2 GeV) = (10.1 ± 1.3 ± 1.3) MeV , (66) 
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where the last error is our estimate of the systematics. One can combine this result with the one for 
the light quark mass ratios from ChPT @: 

r CA _ £bt = 0.553 ± 0.043 , r° A = , 2m ° . = 24.4 ± 1.5 . (67) 

m d ' (m d + m u ) 

Therefore, one can deduce the running masses at 2 GeV: 

m„(2) = (3.6 ±0.6) MeV , m d (2) = (6.5 ± 1.2) MeV , m a (2) = (123.2 ±23.2) MeV . (68) 

Alternatively, we can use the relation between the invariant mass m q and running mass m g (2) to order 
in order to get: 

rh q = (1.14 ±0.05) m q (2) , (69) 
for A3 = (375 ± 50) MeV. Therefore, one can deduce the invariant masses: 

m u = (4.1 ± 0.7) MeV , rh d = (7.4 ± 1.4) MeV , m s = (140.4 ± 26.4) MeV . (70) 

8 Light quark mass from the scalar sum rules 

As can be seen from Eq. ([36|), one can also (in principle) use the isovector-scalar sum rule for 
extracting the quark mass-differences (m d — m u ) and (m s — m u ), and the isoscalar-scalar sum rules 
for extracting the sum (m d + m u ). 

8.1 The scalar ud channel 

In the isovector channel, the analysis relies heavily on the less controlled nature of the ao(980) || 
|55[ |5(^ , [36[ which has been speculated to be a four-quark state |8f|. However, it appears that its qq 
nature is favoured by the present data [86), and further tests are needed for confirming its real qq 
assignemcnt. 

In the 1 = channel, the situation of the ir-ir continuum is much more involved due to the possible 
gluonium nature of the low mass and wide a meson |87], |88|, ^6), which couples strongly to rnr and 
then can be missed in the quenched lattice calculation of scalar gluonia states. 

Assuming that these previous states are quarkonia states, bounds on the quark mass-difference and 
sum of quark masses have been derived in (3^, [l8| [l9| , while an estimate of the sum of the quark masses 
has been recently derived in ]89|] . However, in view of the hadronic uncertainties, we expect that the 
results from the pseudoscalar channels are much more reliable than the ones obtained from the scalar 
channel. Instead, we think that it is more useful to use these sum rules the other way around. Using 
the values of the quark masses from the pseudoscalar sum rules and their ratio from ChPT, one can 
extract their decay constants which are useful for testing the qq nature of the scalar resonances [[| |8q | 
(we shall come back to this point in the next section) . The agreement of the values of the quark masses 
from the isovector scalar channel with the ones from the pseudoscalar channel can be interpreted as 
a strong indication for the qq nature of the ao(980). In the isoscalar channel, the value of the sum of 
light quark masses obtained recently in [ j89[ , though slightly lower, agrees within the errors with the 
one from the pseudoscalar channel. This result supports the maximal quarkonium-gluonium scheme 
for the broad low mass a and narrow /o (980) meson: the narrowness of the fo is due to a destructive 
interference, while the broad nature of the <j is due to a contractive interference allowings its strong 
coupling with tttt. These features are very important for the scalar meson phenomenology, and need 
to be tested further. 

8.2 The scalar us channel 

Here, the analysis is mostly affected by the parametrization of the Kn phase shift data, which strongly 
affects the resulting value of the strange quark mass as can be seen from the different determinations 
given in the Table 0. 
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9 Light quark mass-difference from (M K + — Mk°)qcd 



The mass difference (m^ — m u ) can be related to the QCD part of the kaon mass difference (M K + — 
M k o)qcd from the current algebra relation Q]: 

ca = (™ d - m u ) _ ml (M|p - M 2 K+ ) QCD m\ - to 2 _ 05U0- 3 fr 2 - 1) (711 

2 -(m d + m u )-M| AP K -ml (m u + m d ) 2 ~ ^ ± °-° 5)1 ° ^ 1] ' 

where 2m = to u + m d ; the QCD part of the K + — K° mass-difference comes from the estimate of 
the electromagnetic term using the Dashen theorem including next-to-leading chiral corrections p9[ . 
Using the sum rule prediction of r 3 from the ratio of (m u + m d ) in Eq. ( |66[ ) with the average value of 
m s in Table ^ or the ChPT ratio given in the previous section, one can deduce to order a 3 : 

(m d - Mu) (2 GeV) = (2.8 ± 0.6) MeV. (72) 

Analogous result has been obtained from the heavy-light meson mass-differences |X)| . We shall come 
back to the values of these masses at the end of this chapter. 

10 The strange quark mass from e + e~ and r decays 
10.1 e + e — > I = hadrons data and the 0-meson channel 

Its extraction from the vector channel has been done in [j| and more recently in fi^j , while its 
estimate from an improved Gell-Mann-Okubo mass formula, including the quadratic mass corrections, 
has been done in |)2[ |9^, ||. More recently, the vector channel has been reanalysed in |52| using a 
r-like inclusive decay sum rule in a modern version of the Das-Mathur-Okubo (DMO) sum rule (9^] 
discussed in previous chapter. The analysis in this vector channel is interesting as we have complete 
data from e + e" in this channel, which is not the case of (pseudo) scalar channels where some theoretical 
inputs related to the realization of chiral symmetry have to be used in the parametrization of spectral 
function. One can combine the e + e~ — ► J = 0, 1 hadrons and the rotated recent AS = component 
of the r-decay data in order to extract m s . Unlike previous sum rules, one has the advantage to have 
a complete measurement of the spectral function in the region covered by the analysis. We shall work 
with: 

3|K.| 2 c f M * , A s\ a f 2s\ s 

and the S'C/(3)-breaking combinations p2] |: 

Ai^ = i? Tj i — R T ,(j>, Aio = -R T ,i — 3i? r ,o , (73) 

which vanish in the SU(3) symmetry limit; Aio involves the difference of the isoscalar (R T ,o) an d 
isovector (R T ,i) sum rules a la DMO. The PT series converges quite well at the optimization scale of 
about 1.6 GeV Q. E.g, normalized to m 2 s , one has: 

A 10 ~ -12^|{l + ya s + 30.4a 2 + (173.4 ±109.2)a 3 } 

_LQ«^ Qfi 2 ( m sss - m d dd) 

+ov——r — 3ba„ — . (74) 

The different combinations Ai^ and Aio have the advantage to be free (to leading order) from flavour 
blind combinations like the tachyonic gluon mass and instanton contributions. We have checked using 
the result in |^7| that, to non-leading in to 2 , the tachyonic gluon contribution is also negligible. It has 
been argued in |)5| that Aio can be affected by large SU(2) breakings. This claim has been tested 
using some other sum rules not affected by these terms |5^] but has not been confirmed. The average 
from different combinations is given in Table § An upper bound deduced from the positivity of Rr,<j> 
is also given in Table § 
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10.2 Tau decays 



Like in the case of e + e _ , one can use tau decays for extracting the value of m s . However, data from 
t decays are more accurate than the one from e + e~ . A suitable combination of sum rule sensitive to 
leading order to the SU(3) breaking parameter is needed. It is easy to construct a such combination 
which is very similar to the one for e + e~. One can work with the DMO-like sum rule involving the 
difference between the AS = and AS = —1 processes pq-JrH: 



r>kl Tiki 

5R kl - T - v+A T - s 



|Kd| 2 |K S | 2 
where the moments are defined as: 



^^Ep-e}. (75) 



D>2 



M: 



II 



M*J \M? J ds 
with R®° = R T is the usual r-hadronic width. The QCD expression reads: 



to 2 



(2) 2 (m s ss - m d dd) (4) 



6R» * 24S EW i gAg' - 2 ^-'- M r~' ±% > ■ (77) 



where A^ are perturbative coefficients known to order a 2 : 



^kl — 



i{3Aif| i+r + Aif| L }, (78) 



where the indices T and L refer to the tranverse and longitudinal parts of the spectral functions. 
For D = 2, the L piece converges quite badly while the L + T converge quite well such that the 
combination has can still an acceptable convergence. For the lowest moments, one has: 



A[, 2) = 0.973 + 0.481 + 0.372 + 0.337 
A$ = 1.039 + 0.558 + 0.482 + 0.477 

,(2) 



A { 2 Z > = 1. 115 + 0.643 + 0.608 + 0.647 + ... (79) 

The authors advocate that though the convergence is quite bad, the behaviour of the series is typical 
for an asymptotic series close to their point of minimum snesistivity. Therefore, the mathematical 
procedure for doing a reasonable estimate of the series is to truncate the expansion where the terms 
reach their minimum value. However, the estimate of the errors is still arbitrary. The authors assume 
that the error is given by the last term of the series. The result of the analysis is given in Table ||. The 
different numbers given in the table reflects the difference of methods used to get m s but the results 
are consistent each others within the errors. Like in the case of the e + e~ DMO-like sum rule, the 
combination used here is not affected to leading order by flavour blind contribution like the tachyonic 
gluon and instanton contribution. We have checked that the contribution of the tachyonic gluon to 
order m 2 a s A 2 /M 2 gives a tiny correction and does not affect the estimate done without the inclusion 
of this term. 



10.3 Summary for the estimate of light quark masses 

Here, we summarize the results from the previous analysis: 

• The sum (rn u +rfid) of the running up and down quark masses from the pion sum rules is given 
in Eq. ( p4| ) , while the one of the strange quark mass from the kaon channel is given in Eq. (|5^) . 
Their values lead to the pseudoscalar sum rules prediction for the mass ratio in Eq. ( |57|) which 
agrees nicely with the ChPT mass ratio. 
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• The sum (m u + m<j) of the running up and down quark masses averaged from the pseudoscalar 
sum rule and from a direct extraction of the chiral condensate (uu) obtained from a global fit of 
the nucleon, B* — B mass-splitting and the vector part of the D* — > K*lv form factor is given 
in Eq. @ and reads for A 3 = (375 ± 50) MeV: 

(fn u + m d )(2 GeV) = (10.1 ± 1.8) MeV , (80) 

implying with the help of the ChPT mass ratio m u /md, the value: 

m u (2 GeV) = (3.6 ± 0.6) MeV , m d {2 GeV) = (6.5 ± 1.2) MeV , (81) 

which leads to the invariant mass in Eq. (|7C|): 

rh u = (4.1 ± 0.7) MeV , m d = (7.4 ± 1.4) MeV , (82) 

• We have combined the result in Eq. j6^ ) with the sum rule prediction for m s /(m u + nid) in 
order to deduce the quark mass-difference (mj — m u ) from the QCD part of the K° — K + 
mass- difference. We obtain the result in Eq. ([72]): 

(m d ~ m u ){2 GeV) = (2.8 ± 0.6) MeV . (83) 

This result indeed agrees with the one taking the difference of the mass given previously. The 
fact that (m u + m d ) 7^ {rrid — m u ) disfavours the possibility to have m u = 0. 

• We give in Table || the different sum rules determinations of ra s . The results from the pion 
SR and come from the determination of (m„ + ma) to which we have added the ChPT 
contraint on m s /(m u + md). One can sec from this table that different determinations are in 
good agreement each others. Doing an average of these different results, we obtain: 

m s (2 GeV) = (117.4 ± 23.4) MeV => m s = (133.8 ± 27.3) MeV . (84) 

Aware on the possible correlations between these estimates, we have estimated the error as 
an arithmetic average which is about 10% as generally expected for the systematics of the SR 
approach. 

It is informative to compare the above results with the average of different quenched and unquenched 
lattice values [p6[ : 

m ud (2 GeV) w i(m u + m d )(2 GeV) = (4.5 ± 0.6 ± 0.8) MeV , 

m s {2 GeV) = (110 ± 15 ± 20) MeV , (85) 

where the last error is an estimate of the quenching error. We show in the Table ^ a compilation of the 
lattice unquenched results including comments on the lattice characterisitcs (action, lattice spacing a, 
P). Also shown is the ratio over m s /m u d and quenched (quen) over unquenched (unq) results. 



11 Decay constants of light (pseudo) scalar mesons 
11.1 Pseudoscalar mesons 

Due to the Goldstone nature of the pion and kaon, we have seen that their radial excitations play an 
essential role in the sum rule. This unusual property allows a determination of the radial excitation 
parameters. In the strange quark channels, an update of the results in [^2], ||| ||, |57j gives: 

tk = M K ,f 2 K ,/M K f 2 K ~ 9.5 ± 2.5 ~ rv , (86) 
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Table 4: Simulation details and physical results of unquenched lattice calculations of light quark 
masses from [B6l, where original references are quoted. 



Action a _1 [GeV] 



m s {2) 



_quen 
_unq 



SESAM 98 Wilson 2.3 



MILC 99 Fatlink 1.9 



APE 00 



Wilson 2.6 



CP-PACS 00 MF-Clover a -> 



JLQCD 00 NP-Clover 2.0 



QCDSF + 
UKQCD 00 



NP-Clover 2.0 
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PT 
PT 

NP-PJ 

PT 

PT 
PT 



151(30) {m K ^) 55(12) 1.10(24) 
22(4) 1.08(13) 



113(11) 
125(9) 



112(15) 
108(26) 



88^ 

golf, 

94(2)* 
88(3)* 
109(4)t 
102(6)* 

90(5) 



(m K ) 

K) 



{m K ) 

{m K ) 
( m <p) 



26(2) 1.09(20) 
26(2) 1.25(7) 



{m K ) 
( m <p) 

(m K ) 26(2) 



* From vector WI; * from axial WI. The errors on the ratios m s /m u d and m^ uen /m^ n< ^ are estimates 
based on the original data. 

where rv has been defined previously. The optimal value has been obtained at the LSR scale r 
GeV~ 2 and t c ~ 4.5 - 6.5 GeV 2 as shown in Fig. |. This result implies for tt'(1.3) and K'(1A6): 



f n f ~ (3.3 ± 0.6) MeV , f K , ~ (39.8 ± 7.0) MeV 

It is easy to notice that the result satisfies the relation: 



■fK' M 2 K _ m s 
f n > m 2 m d 

as expected from chiral symmetry arguments. 



(87) 
(88) 



11.2 Scalar mesons 

We expect that the scalar channel is more useful for giving the decay constants of the mesons which 
are not well known rather than predicting the value of the quark masses. Such a programme has been 
initiated in J55|, |5(| |^] . Since then, the estimate of the decay constants has not mainly changed. The 
analysis is shown in Figs. || and ^. Recent estimate gives [52] ]: 



f ao = (1.6 ± 0.15 ± 0.35 ± 0.25) MeV 



f K * ~ (46.3 ± 2.5 ± 5 ± 5) MeV 



(89) 



where the errors are due respectively to the choice of t c from 4.5 to 8 GeV 2 , the value of the quark 
mass-difference obtained previously and the one of A3. The decay constants are normalized as: 



(0\d^(x)\a ) =v / 2/ a M 2 , 



(90) 
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Figure 4: LSR analysis of the ratio r n = M*,f%,/M%f%. For a given value r n — 9.5, we show value of 
(rfid + Su)(2) for two values of the QCD continuum t c . 




corresponding to f n — 92.4 MeV. We have used the experimental masses 0.98 GeV and 1.43 GeV in 
our analysis Q. It is also interesting to notice that the ratio of the decay constants are: 

^29^-^40, (91) 
fa m d ~ m u 

as naively expected. We are aware that the values of these decay constants might have been overesti- 
mated due to the eventual proliferations of nearby radial excitations. Therefore, it will be interesting 
to have a direct measurement of these decay constants for testing these predictions. The values of 
these decay constants will be given like other meson decay constants in the table of the next chapter. 



12 Flavour breakings of the quark condensates 



B The masses of the ao and K£ are also nicely reproduced by the ratio of moments |98j 
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12.1 SU(3) corrections to kaon PC AC 

Let's remind that the (pseudo)scalar two-point function obeys the twice subtracted dispersion relation: 

dt 



*( 5 )(^) = *(5)(0) + g 2 *' (5) (0) + ^ 



t 2 {t 



ie) 



Im*( 5 )(*) 



(92) 



The deviation from kaon PCAC has been firstly studied in |97| using the once subtracted peusdoscalar 
sum rule of the quantity: 



*(5)(9 2 )-*(5)(0) 



(93) 



sensitive to the value of the value of the correlator at q 2 = 10 . The Ward identity obeyed by the 
(pseudo)scalar two-point function leads to the low-energy theorem: 



*(5)(0) = -(mi ± mj)(tpiipi ± ipjipj 



(94) 
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in terms of the normal ordered condensates. However, as emphasized in different papers ]34| . 
\I/(5)(0) contains a perturbative piece which cancels the mass singularities appearing in the OPE 
evaluation of ^(5) (<? 2 ) • This leads to the fact that the quark condensate entering in Eq. H are defined 
as a non-normal ordered condensate, which has a slight dependence on the scale and renormalization 
scheme. This mass correction effect is only quantitatively relevant for the us channel but not for 
the ud one. To order a^ s for the perturbative term and to leading order for the condensates, the 
(pseudo)scalar sum rule for the us channel reads, by neglecting the up quark mass: 

^exp(-frr)-Im* (5 )(f) =s *( 5 )(0) 

t 7T 



+ {m u ± m s ) 2 — t 



-3.15m„r 




(1 - Po) 



1 + 6.82 I — 

7T 



+ 58.55 — 



537.6 — 



(uu) ± (u s) 



10 This sum rule has also been used in 19 9 
which has been checked on the lattice [110 



100|for estimating the U(1)a topological suceptibility and its slope, and 
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+ ^ 2t9 ) (^)^Vs("«) 2 T 3 j, (95) 

where we have neglected the SU (3) breaking for the four-quark condensates. This assumption does 
not however affects the analysis due to the small contribution of this operator at the optimization 
scale. The analysis is shown in Fig. Examining the different curves, on can notice that they deviate 




Figure 7: LSR analysis of the subtraction constant ^5(0). We use m s (2) = 117.4 MeV, tk — 9.5 
and t c = 6 GeV 2 . The curves correspond to different truncations of the PT series: to 0(a s ): dotted- 
dashed; to 0(o? s ): dashed; to 0(a 3 ): continuous. 



notably from the kaon PCAC prediction: 



* 5 (0) ~ 2M 2 K f 2 K , 



(96) 



threfore confirming the early findings in |9j. The LSR indicates a slight stability point at r w (0.50 ~ 
0.75) GeV" 2 , where: 

* 5 (0)~(0.5±0.2)2M2/2 . (97) 
However at this scale, PT series has a bad convergence: 

Pert = Parton x jl + 2.17a s + 5.93a 2 + 17.34a 3 j ~ Parton x |l + 0.86 + 0.92 + 1.06} , (98) 

which might not be worrysome if one considers that asymptotic series close to its point of minimum 
sensitivity can be truncated when its reaches the extremum value and add the last term as truncation 
error |5 This convergence might a priori be improved if one works with the combination of sum 
rules which is less sensitive to the high-energy behaviour of the spectral function (and then to the 
perturbative contribution) than the former |57|, |35], ||, [521. The modified sum rule reads || Q: 

f°° dt 1 3 

/ — exp (-tr) (1 - tr) ~Im%)(i) ~ * (5) (0) + {m u ±m s f —^t' 1 x 

Jo l 7T 07T 



J(« S G 2 ) 



18.3 

8^ 2 



242.2 



— \ 2 
a, 



5.15to 2 t 



[uu) T ^{ss) 



1 + 5.0 
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^ + ^(2T9)(^)7r 3 pa s (^)V 



(99) 



A similar argument has been used for the extraction of the strange quark mass from r-decay data discussed in the 
previous section, where the QCD series also has a quite bad behaviour. 

12 Notice that we have not yet introduced the QCD continuum into the LHS of the sum rule. 
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The analysis also leads to a similar result. The LSR has been also studied recently in [ 102 1 , by 
including threshold effects and higher mass resonances, which enlarge the region of stability in the 
LSR variable. Within the previous hadronic parametrization, one obtains: 



* 5 (0) ~ (0.56 ± 0.04 ± 0.15)2^^/1 , 



(100) 



where we have added the error due to our estimate of the truncation of the QCD PT series as deduced 
from Fig. [?]. An alternative estimate is the uses of FESR [ p2[ . Parametrizing the subtraction constant 
as: 

* 5 (0)« = 2M K f K (l - S K ) , (101) 

one has the sum rule [p2| : 



6k 



m 2 .t r 



167T 2 f%M% 



( M K 



\M K , 



which gives, after using the correlated values of the input parameters |5^, |^, |52) : 



S K = 0.34 



+0.23 
-0.17 ' 



leading to: 

* 5 (0)~ (0.66±0.20)2M|/|- , 
confirming the large violation of kaon PCAC obtained from LSR. 



(102) 

(103) 
(104) 



12.2 Subtraction constant from the scalar sum rule 

One can do a similar analysis for the scalar channel. The analysis from LSR is shown in Fig. 



One 




Figure 8: LSR analysis of the subtraction constant ^(0). We use m s (2) = 117.4 MeV, Jk* 



4G 



MeV and t c = 6.5 GeV . The curves correspond to different truncations of the PT series: to 0(a s ): 
dotted-dashed; to 0{a 2 s ): dashed; to O(a^): continuous. 



can also see that there is a slight stability for r w (0.50 ~ 0.75) GeV 2 , which gives: 

*(0) m -10~ 3 GeV 4 , (105) 
in agreement with previous results U H, HI 13- In P! > 

using LSR, a similar result but from a 
larger range of LSR stability, has been obtained within an Omnes representation for relating the scalar 
form factor to the Kir phase shift data: 



*(0) ~ -(1.06 ± 0.21 ± 0.20)10~ 3 GeV 4 



(106) 
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where the last term is our estimate of the error due to the truncation of the QCD series, the sum rule 
scale One can use an alternative approach by working with FESR: 

= ™\»f K , - j^rn% jl + I a s + 0(a*)} , (107) 

which gives (52]: 

*(0)" = - 10~ 4 GeV 4 . (108) 

12.3 (ss)/(uu) from the (pseudo) scalar sum rules 

We take the arithmetic average of the previous determinations for our final estimate: 

4-5(0) ~ (0.57 ± 0.19)2M£/£ , *(0) ~ -(0.92 ± 0.35)10~ 3 GeV 4 , (109) 

Taking the ratio of the scalar over the pseudoscalar subtraction constants expressed in terms of the 
normal- ordered condensates, one can deduce: 

(ss)/{uu) = 0.57 ±0.12. (110) 

12.4 (ss)/(uu) from the B s meson 

One can also extract the flavour breakings of the condensates from a sum rule analysis of the B s and 
B* masses, which are senstive to the chiral condensate as it enters like mi,(ss) in the OPE of the 
heavy light meson (see next section). The masses of the mesons are found to decrease linearly with 
the value of the chiral condensate. Using the observed value of the B s meson mass Mb s = 5.375 GeV, 
one can deduce from Fig 3 of pOf : 

(ss)/(uu) -0.75 ±0.08 , (111) 

where the error is the expected typical sum rule estimate. The effect of the strange quark mass is 
less important than this one here, such that the result given in [|so| remains valid though obtained 
with slightly different values of m s and A5. This estimate is expected to be more reliable than the 
one from the (pseudo)scalar light mesons, which are affected by the bad convergence of the PT QCD 
series. Using this value of ratio of the condensates in the curve of the B* mass, one can also predict: 

B* s = 5.64 GeV , (112) 

which can be tested in B factories. 



12.5 Final sum rule estimate of (ss)/{uu) 



Using the previous results, one can deduce that the sum rules from the light (pseudo)scalar and from 
the B s meson predict for the normal ordered condensate ratio: 

(ss)/(uu) ~ 0.66 ±0.10 , (113) 

on the large flavour breaking of the chiral condensate. This 



confirming earlier findings [ p7[ [3| p5| , 

number comes from the arithmetic average of the two values in Eqs. (110) and ( 



111 



). If one instead 

works with the non-normal ordered condensate, one should add to the expression in Eq. ( |94| ) a small 
perturbative piece first obtained by Becchi et al (35) (see also B |{|3, |4p[|): 



(ss) 



MS 



3 2 
2^2 7 



1 



53 
24 



mi 



(114) 



24 



This leads to the ratio of the non-normal ordered condensates: 



(ss 



)/{uu)\ m = 0.75 ±0.12. (115) 



The prev ious estimates are in good agreement with the ones from chiral perturbation theory (see 



also [ 103 1 ) . They are also in fair agreement with the one from the baryonic sum rules [fT2|-|78(|, though 
we expect that the result from the latter is less accurate due to the complexity of the analysis in 
this channel (choice of the interpolating operators, eventual large effects of the continuum due to the 
nearby Roper resonances,...). 

12.6 SU(2) breaking of the quark condensate 

The SU(2) breaking of the quark condensate has been studied for the first time in ^] and in [[32| [|. 
Using similar approaches, the estimate is ||, |5f|: 

(dd)/(uu) ~ 1 - 9 x 10" 3 . (116) 

The previous estimate is in good agreement with the one from FESR []62[ . 

13 Heavy quark masses 

In the previous part of this book, we have already discussed the different definitions of the heavy 
quark masses and given their values. Contrary to the light quark masses, the definition of pole quark 
masses p 2 — Mfj can (in principle) be introduced perturbatively for heavy quarks Jll], |l2| [l3| similarly 
to the one of the electron as here the quark mass is much heavier than the QCD scale A such the 
perturbative approach makes sense. However, a complication arises due to the resummation of the 
QCD series [^3 such that the pole mass definition has an intrinsic ambiguity, which can be an obstacle 
for its improved accurate determination though the effect is relatively small. Alternative definitions 
free from such ambiguities have been proposed in the literature . In this section, we shall discuss 

the determinations of the perturbative running quark masses which do not have such problems. 

13.1 The quarkonia channel 

Charmonium and bottomium are the standard channels for extracting the charm and bottom quark 
masses. Most of the sum rule analysis are based on the Q 2 = moments (MOM) originally introduced 
by SVZ for the study of the charmonium systems: 



7^T ilmn « . ( 117 ) 

4m 2 t n+i 7T 



but convenient for the bottomium systems due to a much better convergence of the OPE. In [981, the 
Q 2 7^ moments have been introduced for improving the convergence of the QCD series: 



1 



dt 1 



4m 2 (* + QIY+ 1 7T 



The spectral function can be related to the e + e — > QQ total cross-section via the optical theorem: 

1 a(e+e- -» QQ) 

Imn t + ze = - \ _ ^ . 119 

YI-kQ^q a(e+e — > (i+fi ) 

Qq is the heavy quark charge in units of e. The contribution to the spectral function is as usual 
saturated by the lowest few resonances plus the QCD continuum above the threshold t c : 

ImR Q (t) = ^-± r ^2T i M i 6(t-M?) + 6(t-t c )ImIl% CD (t), (120) 
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where 1^ is the electronic width of the resonances with the value given in PDG pOfl . Retaining the 
observed resonances, the value of ^Jt~ c fixed from stability analysis is about (11 ~ 12) GeV for the T- 
and about 5 GeV for the -families. However, the result will be practically independent from this 
choice of t c due to the almost complete dominance of the lowest ground state to the spectral function 



at the stability point. An alternative approach used in |13j, 104 1 is the LSR: 

f°° 1 

C(t) = / dt exp-* T -Imn(t) . (121) 
This sum rule is particularly convenient for the analysis of the charmonium systems as the correspond- 



ing OPE converges faster than the moment sum rules. It has been noticed in [104] that the ratios of 
sum rules (and their finite energy sum rule (FESR) variants) are more appropriate for the estimate of 
the quark mass as these ratios equate directly the mass squared of ground state to that of the quark: 

M {n) d 
T^ = J^ and K T = --\agC, (122) 

They also eliminate, to leading order, some artefact dependence due to the sum rules (exponential 
weight factor or number of derivatives) and some other systematic errors appearing in each individ- 
ual moments. For the perturbative part, we shall use (without expanding in 1/M) the Schwinger 
extrapolation formula to two-loops: 



Imn^(t) ~ ±-v Q (^-pA {l + ^ S /M} , (123) 



where: 



(3 + vq) fn 3 



V Q = y/l-UQ/t , f(v Q ) = JL - ^ - £ ) 



are respectively the quark velocity and the Schwinger function 1 105 . We express this spectral function 
in terms of the running mass by using the two-loops relation given in previous chapter and including 
the a s log(t/Mg)-term appearing for off-shell quark. We shall add to this perturbative expression the 
lowest dimension (a s G 2 ) non-perturbative effect (it is known as explained in the previous part of the 
book that, for a heavy-heavy quark correlator, the heavy quark condensate contribution is already 
absorbed into the gluon one) which among the available higher dimension condensate-terms can only 
give a non-negligible contribution. The gluon condensate contribution to the moments MS n ' and so 
to lZ n can be copied from the original work of SVZ pT) and reads: 



,(») _ K/I {n) (n + 3)! Air (a s G 2 
(n-l)!(2n + 5) 9 f^y 



(n) 

where M. pert is the lowest perturbative expression of the moments. The one to the Laplace ratio 1Z T 
can be also copied from the original work of Bertlmann |106| , which has been expanded recently in 
l/M Q by |107| . It reads: 

K^ {i M^ {a ,Giy-( i + ±-— 2 ), (126) 

where u> = AMqT. The results of the analysis from the ratios of moments and Laplace sum rules give 
the values of the running masses to order a s [^] : 

m c (m c ) = (1.23 ± 0.03 ± 0.03) GeV , m b {m b ) = (4.23 ± 0.04 ± 0.02) GeV , (127) 



13 The inclusion of the correction is under study. 
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where the errors are respectively due to a s {M z ) — 0.118±0.006 and (a B G 2 ) = (0.06±0.03) GeV 4 used 
in the original work. These running masses can be converted into the pole masses at this order. Non- 
relativistic versions of these sum rules (NRSR) introduced by | 10£ | have also been used in 104 for 



determining the b quark mass. These NRSR approaches have been improved by the inclusion of higher 
order QCD corrections and resummation of the Coulomb corrections from ladder gluonic exchanges. 
Some recent different determinations are given in Tables | and | 



Table 5: QSSR direct determinations of m c (rn c ) in MS'-scheme and of the pole mass M c from J/\f'- 
family, e + e— data and Z?-meson and comparisons with lattice results. Determinations from some 
other sources are quoted in PDG ]2C[ ]. The results are given in units of GeV. The estimated error in 
the SR average comes from an arithmetic average of the different errors. The average for the pole 
masses is given at NLO. The one of the running masses is almost unchanged from NLO to NNLO 
determinations. <= means that perturbative relation between the different mass definitions have 
been used to get the quoted values. 



Sources 



m c (m c ) 



Comments 



Authors 



J/\I>-family 

MOM and LSR at NLO (1.27 ± 0.02) 



Ratio of LSR at NLO 
NRSR at NLO 
SR at NLO 
NRSR at NNLO 



(1.23 ±0.04) 
(1.23 ±0.04) 
(1.22 ±0.06) 
(1.23 ±0.09) 



(1.45 ±0.05) 
(1.42 ±0.03) 
(1.45 ±0.04) 
(1.46 ±0.04) 
(1.70 ±0.13)* 



(1.26 ±0.02) 



SN89 pi 
SN94 [104 1 
SN94 jicjl 
DGP9 4 1107 ] 
EJ01 Jll2t 



e+e data 

FESR at NLO 
MOM at NNLO 
NLO 



(1.37 ±0.09) 
(1.30 ±0.03) 
(1.04 ±0.04) 



1.33- 1.4 



PS01 hit. 
KS01 ITu 



M01 [115] 



D meson 

Ratio of LSR at NNLO (1.1 ± 0.04) (1.47 ± 0.04) 



SN01 flll 



SR Average 



(1.23 ±0.05) (1.43 ±0.04) 



Quenched Lattice 



(1.33 ±0.08) 
(1.20 ±0.23) 
(1.26 ±0.13) 



FNAL98 Hi 16 1 
NRQCD99 



APE01 [118 



117 1 



* Not included in the average. 



27 



13.2 The heavy-light D and B meson channels 

Heavy quark masses can also be extracted from the heavy-light quark channels because the corre- 



sponding correlators are sensitive to leading order to the values of these masses |3j, |80|, 104, 109, 110|. 
Again, we shall be concerned here with the LSR C(t) and the ratio TZ(t). The latter sum rule, or its 
slight modification, is useful, as it is equal to the resonance mass squared, in the simple duality ansatz 
parametrization of the spectral function: 



1 



In# 5 (t) ~ f D M%6{t - M%) + "QCD continuum" 6(t-t c ), 



(128) 



where fr> is the decay constant analogue to / w = 130.56 MeV. The QCD side of the sum rule reads: 

1 



£qcd(t) 



A0 



dt 



-Mr 



8tt 2 



R2s 



C 4 (Q 4 ) + C 6 (0 6 )n 



(129) 



where R2s is the new a 2 -term obtained semi-analytically in [119| and is available as a Mathematica 
package program Rvs.m. Neglecting md, the other terms are: 

X = Mg/t, 

f(x) = ^+2Li 2 (x)+logxlog(l-x)-^log(l/x-l) 



C 4 (0 4 ) 

c 6 (o 6 ) 



— log(l — x) + a;log(l/x — 1) — (x/(l — x)) log a;, 
-M Q (dd) + (a s G 2 ) /12tt 

3tt\ / Mir Mir 2 \ _ „ 



(130) 



The previous sum rules can be expressed in terms of the running mass toq [y) |^[ From this expression, 
one can easily deduce the expression of the ratio H(t), where the unkown decay constant disappears, 
and from which we obtain the running quark masses: 



m c (m c ) = (1.10 ±0.04) GeV 



(131) 



The analysis is shown in Fig. [)], where a simulataneous fit of the decay constant from C and of m c (m c ) 
from 1Z is shown |^[ 



Our optimal results correspond to the case where both stability in r and in t c are reached. However, 
for a more conservative estimate of the errors we allow deviations from the stability points, and we 
take: 

t c ^ (6 - 9.5) GeV 2 , t ~ (1.2 ± 0.2) GeV -2 , (132) 

and where the lowest value of t c corresponds to the beginning of the r-stability region. Values outside 
the above ranges are not consistent with the stability criteria. One can inspect that the dominant 

14 Notice that we have adopted here the lattice normalization for avoiding confusion. We shall discuss its determination 
in the next chapter. 

15 It is clear that, for the non-perturbative terms which are known to leading order of perturbation theory, one can 
use either the running or the pole mass. However, we shall see that this distinction does not affect notably the present 
result. 

16 We shall discuss the decay constant in the next section. 
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t[GeV z ] S 6 (m„) [GeV] 



Figure 9: Laplace sum rule analysis of fn and m c (m c ). 



non-perturbative contribution is due to the dimension- four M c (dd) light quark condensate, and test 
that the OPE is not broken by high-dimension condensates at the optimization scale. However, 
the perturbative radiative corrections converge slowly, as the value of fp increases by 12% after 
the inclusion of the a s correction and the sum of the lowest order plus a s -correction increases by 
21 % after the inclusion of the a 2 term, indicating that the total amount of corrections of 21% is 
still a reasonnable correction despite the slow convergence of the perturbative series, which might be 
improved using a resummed series. However, as the radiative corrections are both positive, we expect 
that this slow convergence will not affect in a sensible way the final estimate. A similar analysis is 
done for the pole mass. The discussions presented previously apply also here, including the one of the 
radiative corrections. We quote the final result: 

M c = (1.46 ± 0.04) GeV , (133) 

where the error is slightly smaller here due to the absence of the subtraction scale uncertainties. One 
can cross-check that the two values of fh c (m c ) and M c give the ratio: 

M c /fh c (m c ) ~ 1.33 , (134) 

which satisfies quite well the three-loop perturbative relation M c /fh c (rn c ) = 1.33. This could be a non- 
trivial result if one has in mind that the quark pole mass definition can be affected by non-perturbative 
corrections not present in the standard SVZ-OPE. In particular, it may signal that l/q 2 correction of 
the type discussed in [^(| |2g}] , if present, will only affect weakly the standard SVZ-phenomenology 
as observed explicitly in the light quark, gluonia and hybrid channels p7ft . Using analogous analysis 
for the B meson, we obtain at the optimization scale r = 0.375 GeV" 2 and t c — 38 GeV 2 : 

fh b {m b ) = (4.05 ± 0.06) GeV , (135) 

while using the pole mass as a free parameter, we get: 

M b = (4.69 ± 0.06) GeV , (136) 
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One can again cross-check that the two values of mb(mb) and M b lead to 

Mb/fhb(mb) = 1.16 , 



(137) 



to be compared with 1.15 from the three-loop perturbative relation between Mb and nib, an d might in- 
directly indicate the smallness of the l/q 2 correction if any. One can immediately notice the agreement 
of the results from quarkonia and heavy-light quark channels. Comparisons with other determinations 
are given m Tables § and |. 

Table 6: The same as in Table pi but for the 6-quark. 



Sources 



m b (m b ) 



Mb 



Comments 



Authors 



Y-family 

MOM and LSR at NLO 
Ratio of LSR at NLO 
NRSR at NLO 
FESR at NLO 

NRSR at NNLO 
MOM at NNLO 
NR at NNNLO 
NR at NNNLO 
NR at NNLO 
NR at NNLO 
MOM at NNNLO 

B and B* mesons 

Ratio of LSR at NLO 
Ratio of LSR at NNLO 



(4.24 ±0.05) 
(4.23 ±0.04) 
(4.29 ±0.04) 
(4.22 ±0.05) 
(4.14 ±0.04) 
(4.20 ±0.10) 
(4.19 ±0.06) 
(4.45 ±0.04) 
(4.21 ±0.09) 
(4.25 ±0.08) 
(4.20 ±0.06) 
(4.21 ±0.05) 



(4.24 ±0.07) 
(4.05 ±0.06) 



(4.67 ±0.10) 
(4.62 ±0.02) 
(4.69 ±0.03) 
(4.67 ±0.05) 
(4.75 ±0.04) 



(4.63 ±0.08) 
(4.69 ±0.06) 



m b (- 



= (4.23 ±0.05) SN89 (T| 



Residual mass 
IS mass 



SN94 104 



SN94 1 104 



SN95 (104 



KPP98 [|120| 



PP99, MY99 |21| 

JP99 m 



PY00, LS00 (123 1 
P01 p| 
BS99 J125] 
H00 p6j 
KS01 TU4 



_B-meson only 



SN94 
SN01 



1 104 
[110 



SR Average 



(4.24 ± 0.06) (4.66 ± 0.06) =>■ m b {M z ) = (2.83 ± 0.04) 



Average LEP 

3-jets at Mz 

Unquenched Lattice 



(4.23 ± 0.94) 



(4.23 ±0.09) 



m b (Mz) = (2.82 ± 0.63) LEP |127] 



APE00 1 128 



Summary for the heavy quark masses and consequences 

From Tables | and |, we conclude that the running c and b quark masses to order a 2 s from the different 
sum rules analysis are likely to be: 

m c (m c ) = (1.23 ± 0.05) GeV , 



m b {m b ) = (4.24 ± 0.06) GeV , 



(138) 
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where the estimated errors come from the arithmetical average of different errors. Wc have not tried to 
minimize the errors from weighted average as the correlations between these different determinations 
are not clear at all. However, as one can see in the tables, the quoted errors are typical for each 
individual determinations. These results are consistent with other determinations given in |20| and 
in particular with LEP average from three-jet events and lattice values reported in the tables. Using 
the previous relation between the short distance perturbative pole and running masses, one obtains, 
to order a s : 



M. 



PT2 



(1.41 ±0.06) GeV , 



and to order ai: 



M. 



PT3 



= (1.64 ±0.07) GeV , 



M b PT2 



M b PT3 



(4.63 ±0.07) GeV 



(4.88 ±0.07) GeV 



(139) 



(140) 



which are consistent with the average values to order a s quoted in the tables and in 104]. However, 
one should notice the large effects due to radiative corrections which can reflect the uncertainties in 
the pole mass definition. From the previous values of the running masses, one can also deduce the 
values of the RG invariant masses to order al : 



(1.21 ±0.07) GeV 



m b = (6.9 ±0.2) GeV 



(141) 



We have used A 4 = 325 ± 40 MeV and A 5 = 225 ± 30 MeV. Taking into account threshold effects and 
using matching conditions, we can also evaluate the running masses at the scale 2 GeV and obtains: 



m c (2) = (1.23 ±0.05) GeV , 



m b (2) = (5.78 ±0.26) GeV , 



(142) 



Combining the values of mb and m s obtained in the previous section, one can deduce the scale 
independent mass ratio: 

I J ! ; 

(143) 



-- =48.8 ±9.8 , 



which is useful for model buildings. 

One can also run the value of mj at the Z-mass, and obtains the value of rrib(Mz) quoted in the table: 



m b (M z 



(2.83 ±0.04) GeV. 



(144) 



This value compares quite well with the ones measured at Mz from three-jet heavy quark production 
at LEP where the average (2.83 ± 0.04) GeV of different measurements ]127[ is given also in the 
table. This is a first indication for the running of mb in favour of the QCD predictions based on the 
renormalization group equation. 



14 The weak leptonic decay constants and fs {s) 

In this section we summarize the different results obtained from the QCD spectral sum rules 
(QSSR) on the leptonic decay constants of the B and D mesons which are useful in the analysis of the 
leptonic decay and on the B-B mixings. Intensive activities have been devoted to this subject during 
the last few years using QSSR and lattice calculations. 
The leptonic constant of the pseudoscalar P = D, B meson is defined as: 

(0\d^\P) = fpMpP , (145) 

where P is the pseudoscalar meson field and fp is the pseudoscalar decay constant which controls the 
P — > Iv leptonic decay width, nor malized as f„ = 130.56 MeV0- The current : 

d^jx)) = (mi + M^tpiji'j^ipj (i = u, d, s; j = c, b) , (146) 

17 This is an extension and an update of the some parts of the reviews given in g. 

18 In this chapter, we adopt this normalization used by the lattice and experimental groups. In the previous sections, 
we have used / x = / w /v2 MeV. 
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is the divergence of the axial current. In the sum rule analysis, we shall be concerned with the 
pseudoscalar two-point correlator: 

* 5 (g 2 ) = i /Vze^<0|Tc^(x)} (c^(0);.) f |0) . (147) 

In the case of the B(ub) meson, the decay width into tv t reads: 

T(B - ru T + B - rv Tl ) = ^f^M B (l - ' Af T 2 / 2 , (148) 

where M T expresses the helicity suppression of the decay rate into light leptons e and \i. This expression 
shows that a good determination of fs will allow a precise extraction of the CKM mixing angle V u b- 
One the other hand, fs and the so-called bag parameter Bb also control the matrix element of the 
AB = 2 B°-B° mixing matrix, which is of a non-perturbative origin, as we shall discuss in another 
chapter. 

However, contrary to the case of the tt and K mesons, the leptonic width of the heavy meson is small 
as the corresponding decay constant vanishes as 1/ J Mq , while the presence of the neutrino in the 
final state renders difficult the reconstruction of the signal and the rejection of background. Moreover, 
the B leptonic rate is Cabibbo suppressed, which makes it unreachable with present measurements. 
(~ |Kib| 2 ), while the D s leptonic rate is Cabbibbo favoured (~ |V CS | 2 ). Recent measurements of fn s 
are given in Fig. [l(], where the quoted average is |129| : 

f Ds ~ (264 ± 37) MeV . (149) 
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Figure 10: Different measurements of fn s compared with theoretical predictions from [ 129 1 . 



14.1 Upper bound on the value of fn 



Wi thin the QSSR framework, the decay constants of the B and D mesons have been firstly estimated 
in 1 131], while the first upper bounds on their values have been derived in [132] and updated in 
the recent review 0]. Indeed, a rigorous upper bound on these couplings can be derived from the 
second-lowest superconvergent moment: 



_ 1 d 2 ^ 5 (q 2 



2! 



(dq 2 f 



(150) 
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where for this low-moment, the OPE well beh aves. Using the positivity of the higher-state contribu- 
tions to the spectral function, one can deduce [132, |93|j : 



M P 
fp< — 



M, 



0.751a s + .. 



Q 



(151) 



where one should not misinterpret the mass-dependence in this expression compared to the one ex- 
pected from heavy-quark symmetry. Applying this result to the D meson, one obtains: 



Id < 2.14/, 



(152) 



which is not dependent to leading order on the value of the charm quarkm mass. Although presumably 
quite weak, this bound, when combined with the recent determination of the SU(3)p breaking effects 
to two loops on the ratio of decay constants [133]: 



implies 



(1.15 ±0.04)/, , 

ID 



f D , < (2.46 ± 0.09)/, ~ (321.2 ± 11.8) MeV 



(153) 



(154) 



which is useful for a comparison with the recent measurement of fo 3 , with the average value given 
previously. One cannot push, however, the uses of the moments to higher n values in this D channel, 
in order to minimize the continuum contribution to the sum rule with the aim to derive an estimate 
of the decay constant from this method, and to derive its "correct" mass dependence, because the 
QCD series will not converge at higher n values. 



14.2 Estimate of the D decay constant fjj 



The decay constant fp> can be extracted from the pseudoscalar Laplace sum rules given in Eq. ( |129| ) 
0. Prior 1987, the different sum rules estimate of the decay constant fp have been inconsistent among 
each others. To our knowledge, the first attempt to understand such disrepancies has been done in 



[134| 135] (see also [ 136 1 ) , where it has been shown, for the first time and a long time before the lattice 
results, that: 

Id « f B « (1.2 ~ 1.5)/, , (155) 



which differs from that expected from the heavy quark symmetry scaling law [ 137 1 : 



fB 



M D f a s {M c ) 
Mp \a s (M b ) 



-1/01 



(156) 



valid in the extreme case where the heavy quark mass is infinite f^j It has also been understood 
that the apparent disagreement among different existing QSSR numerical results in the literature is 
not only due to the choice of the continuum threshold t c [ its effect is (7 ~ 10)% of the result when 
one moves t c from the one at the beginning of sum rule variable to the one where the t c stability 
is reached. ]q] as misleadingly claimed in the literature. Indeed, the main effect is also due to the 
different values of the quark masses used which is shown explicitly in the table of 133 ] . 
In the D channel, the most appropriate sum rule is the (relativisitc) Laplace sum rule, as the OPE 
of the q 2 = moments does not converge for larger value of the number of derivatives n, at which 
the D meson contribution to the spectral integral is optimized. The results from different groups are 



i9p or reviews, see e.g. [130 . 

20 Finite mass corrections to this formula will be dicussed later on. 

21 In some papers in the literature, the value of t c is taken smaller than the previous range. In this case, the t c effect 
is larger than the one given here. 



1>2 



A critical review on the discrepancy between different existing estimates is given in [ 13S ] 
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consistent with each others for a given value of the c-quark mass. For the D meson, the optimal result 
is obtained for: 

6 < t c < 9.5 GeV 2 , r ~ (1.2 ± 0.2) GeV" 2 . (157) 

where the QCD corrections are still reasonnably small. The most recent estimate including a? cor- 
rections from a simultaneous fit of the set either (fD,rn c (m c )) or (fD,M^ ole ) is given in Fig. p[ Th e 
obtained values of the quark masses have been quoted in Table [|. The resulting value of fr> is |llCfl : 

f D ~ (203 ± 23) MeV , (158) 

in ag reement with the recent evaluation (195 ± 20) MeV at order a 2 but using the pole mass as input 
fl39f. 



14.3 Ratio of the decay constants fo s l fo an d fsj f~B 



The 5J7(3) breaking ratios /d 3 //d and /b s //b have been obtained semi-analytically in [ 133 1 . In 
order to have a qualitative understanding of the size of these corrections, we strat from the global 
hadron-quark duality sum rule: 



dw Im*f Q (w) 



dw Im% es (uj) : 

where uj c is the continuum energy defined as: 

t = (E + M Q f = M% + ujM, 
Keeping the leading order terms in a s and in 1/Mq, it leads to 



R P ~p P {l 



1 - 



Mr 



- 6 



where: 



Pp 



_ ( M P \ 



MpJ V 



ms 2 

LO c 

TO, 

H 

Mr 



Mr 



1 - 



Mr 



The value of u> is fixed from stability criteria to be 134, 138, 141, 142 

w c ~ (3.1 ±0.1) GeV. 

The sum rule indicates that the SU(3) breaking corrections are of two types, the one to s /Mq and the 
other m s /uj c . More quantitatively, we work with the Laplace sum rule: 



(159) 
(160) 

(161) 

(162) 
(163) 



L = 

Analogously the Laplace sum rule gives: 
Rp~pp \ 1 + 2(2.2 ±0.2) 



m s 
M~r 



2 (8.2 ± 1.6) 



ms 2 



(164) 



(165) 



where the numerical integration includes a slight Mq dependence. Including m s a s and m 2 a s - 
corrections, the resulting values of the ratio are: 

R D = ^ = 1.15 ± 0.04 , R B = ^ = 1.16 ± 0.05 . (166) 

JD IB 

This result implies: 

fx) g ~ (235 ± 24) MeV , (167) 
which agrees within the errors with the data [ 129 1 and lattice ]140[ averages both quoted in Fig. |l^. 
This feature increases the confidence in the uses of the QSSR method for predicting the not yet 
measured decay constant of the B meson. 
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Table 7: Estimate of /s (s) to order a 2 s and /s»//b to order a s from QSSR and comparison with the 
lattice. 



Sources 


fs [MeV] 


/s (s) / Jb 


fs 3 [MeV] 


Comments 


Authors 


QSSR 

LSR 

HQETSR 


203 ± 23 
210 ± 19 
206 ± 20 


1.16 ±0.04 = 


236 ± 30 
244 ± 21 


M po i e , m b : output 
mb: input 
M po i e : input 


SNOl [110| 
JLOl |l4| 

psoi psft 


SR Average 


207 ±21 =^ 




240 ± 24 






Unq. Lattice 


200 ± 30 


1.16 ±0.04 = 


> 232 ± 35 


average 


LAT01[144| 



14.4 Estimate of the B decay constant f B 

For the estimate of fs, one can either work with the Laplace, the moments or their non-relativistic 
variants because the 6-quark mass is relatively heavy. The optimal result which we shall give here 
comes from the Laplace relativistic sum rules. They corresponds to the conservative range of param- 
eters: 

0.6 < = - Mb < 1-8 GeV , r ~ 0.38 GeV~ 2 , (168) 



which have been used in the previous section for getting the 6-quark mass 
the dominant corrections come from the ( 



As shown in the figure 
quark condensate with the strengh 



given in [ 14£ , 146 

(30 ~ 40)% of the lowest order term in fg, while the higher condensate effects are smaller, which are 
respectively —(20 ~ 30)% , + (5 ~ 8)% for the d — 5 and 6 condensates. This shows, despite the 
large value of the quark condensate contribution, that the OPE is convergent. It has been noticed in 
145, |146[ , that the convergence of the OPE is faster for the relativistic LSR than for the moments, 
such that the most precise result should from the LSR . In both cases the perturbative corrections are 



small. One obtains from the relativistic LSR , the results to order a 2 [ 1 1 1 

f B ^ (203 ± 23) MeV ~ (1.55 ± 0.18)/^ , 
and to order a s (see previous discussion) ]133j : 

~ 1.16 ±0.04 . 



/b. 

fB 



(169) 



(170) 



These values of fg and fp agree quite well with the previous QSSR findin gs in |133| , and [136 



They also agree with the non-relativistic sum rules estimate in the full theory [134], in HQET [147, 1 



and in 141 142 



However, unlike the relativistic sum rules, the HQET sum rule is strongly affected 
by the huge perturbative radiative corrections of about 100%, which is important at the sum rule scale 
of about 1 GeV at which the HQET sum rule optimizes. These results are also in good agreement 
with the lattice average estimate given in Table [t| 
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14.5 Static limit and l/M b -corrections to f B 



As noticed previously, the first result fg ~ fjj in [133], which has been confirmed by recent estimates 
from different approaches, shows a large violation of the scaling law expected from heavy-quark sym- 
metry. This result suggests that finite quark mass corrections are still huge at the D and B meson 
masses. The first attempt to understand this problem analytically is in ]135| in terms of large cor- 
rections of the type E c /M{> if one expresses in this paper the continuum threshold t c in terms of the 
threshold energy E c : 

t c = (E c + M b ) 2 . (171) 

Later on different approaches have been investigated for the estimate of the size of these corrections. 
In the lattice approach, these mass corrections have been estimated from a fit of the obtained value 
of the meson decay constant at finite and infinite (static limit) values of the heavy quark mass and 
by assuming that these corrections are polynomial in 1/Mn up to log. corrections. A similar analysis 
has been done with the sum rule in the full theory |l45| , |l46| , by studying numerically the quark mass 
dependence of the decay constant until the quark mass value [Mq < 15 GeV) until which one may 
expect that the sum rule analysis is still valid. In so doing, we use the parametrization: 



.4 



B 



M h Mi 



L 3 7T 

by including the quadratic mass corrections. The analysis gives 0: 

Ib = (/ba/a/b) - (0.65 ± 0.06) GeV 3 / 2 



(172) 



(173) 



which one can compare with the results from the HQET Laplace sum rule [ 147 1 and [142, |148|] : 

f B ^ (0.35 ± 0.10) GeV 3/2 , (174) 



and from FESR JL4J: 



] B ^ (0.57 ±0.10) GeV 3/2 , 
Taking the average of these three (independent) determinations, one can deduce: 

f B ~ (0.58 ± 0.09) GeV 3/2 , 



(175) 



(176) 



where we have done an arithmetic average of the errors. This result is in good agreement with the 
lattice value given in [14C, 149 using nonperturbative clover fermions. One can translate this result 
into the value of fg in the static limit approximation: 



/ B * Qt ~(1.9±0.3)/ w . 



(177) 



We can also use the previous value of f B tat together with the previous values of fg and fu at the 
"physical" quark masses in order to determine numerically the coefficients A and B of the 1 /Mb and 
1/Af 2 corrections. In so doing, we use the values of the quark "pole" masses given in Tables | and §. 
Then, we obtain from a quadratic fit: 



A re 0.98 GeV and B ps 0.35 GeV 2 



(178) 



while a linear fit gives a large uncertainty: 



A& (0.74-0.91) GeV 



(179) 

23 The numbers given in [145, |l46| correspond to the quark mass Mf, = 4.6 GeV and should be rescaled until the 
meson mass Mg. In the following, we shall also work with fg normalized to be \/2 bigger than in the original papers. 
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One can notice that the fit of these coefficients d epends s t rong ly one the input values of fry and 
fs- Indeed, using some other set of values as in [145, 146, 141], one obtains values about 2 times 
smaller. Therefore, we consider as a conservative value of these coefficients an u ncert ainty of about 
50%: The value of A obtained is comparable with the one from HQET sum rules 150 and 145 , 141] 
of about 0.7 ~ 1.2 GeV. A similar value of A has been also obtained from the recent NR lattice 



calculations with dynamical fcrmions and using a linear fit [149 

A w 0.7 GeV . 



(180) 



One can qualitatively compare this result with the one obtained from the analytic expression of the 
moment in the full theory H, fTM p|: 



l_El 

TT 2 M b 



f 2 ~ 

JB ~ _9 n,r 



1 _2o ! _ 3{n + 2)E c 

3 7T M b 



2 m 



(181) 



Here, one can notice that the size of the l/Mft corrections depends on the number of moments, such 
that their estimate using literally the expression of the moments can be inconclusive. A qualitative 
estimate of these corrections can be done from the semilocal duality sum rule, which has more intuitive 
physical meaning due to its direct connection with the leptonic width and total cross-sec tion through 
the optical theorem. It corresponds to n = —2, and leads to the interpolating formula [148]: 




3 7t 88 Ml 



2 R3 



(182) 



from which, one obtains: 

A w ^(Mg — Mb) ~ 1 GeV, 
3 27 

B » -£ c 2 + -(M B -A/ b ) 2 ~ 0.45 GeV 2 , 
which is in good agreement with the previous numerical estimate. 



(183) 



15 Conclusions 

We have reviewed the determinations of the quark masses and leptonic decay constants of (pseudo)scalar 
mesons which are useful in particle physics phenomenology. The impressive agreements of the QSSR 
results with the data when they are measured and/or with lattice calculations in different channels in- 
dicate the self-consistency of the QSSR approach making it as one of the most powerful semi-analytical 
QCD nonperturbative approach available today. Applications of these results for studying the B-B 
mixings and CP- violation will be discussed in the next chapter. 
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